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Abstract 

Nentrino spin and flavour oscillation in curved spacetime have been studied for the most 
general static spherically symmetric configuration. Having exploited the spherical symmetry 
we have confined ourselves to the equatorial plane in order to determine the spin and flavour 
oscillation frequency in this general set-up. Using the symmetry properties we have derived 
spin oscillation frequency for neutrino moving along a geodesic or in a circular orbit. Starting 
from the expression of neutrino spin oscillation frequency we have shown that even in this 
general context, in high energy limit the spin oscillation frequency for neutrino moving along 
circular orbit vanishes. We have verified previous results along this line by transforming to 
Schwarzschild coordinates under appropriate limit. This finally lends itself to the probability 
of neutrino helicity flip which turns out to be non-zero. While for neutrino flavour oscillation 
we have derived general results for oscillation phase, which subsequently have been applied 
to three different gravity theories. One, of them appears as low-energy approximation to 
string theory, where we have an additional held, namely, dilaton held coupled to Maxwell 
held tensor. This yields a realization of Reissner-Nordstrom solution in string theory at low- 
energy. Next one corresponds to generalization of Schwarzschild solution by introduction 
of quadratic curvature terms of all possible form to the Einstein-Hilbert action. Finally, we 
have also discussed regular black hole solutions. In all these cases the flavour oscillation 
probabilities can be determined for solar neutrinos and thus can be used to put bounds 
on the parameters of these gravity theories. While for spin oscillation probability, we have 
considered two cases, Gauss-Bonnet term added to the Einstein-Hilbert action and the 
f(R) gravity theory. In both these cases we could impose bounds on the parameters which 
are consistent with previous considerations. In a nutshell, in this work we have presented 
both spin and flavour oscillation frequency of neutrino in most general static spherically 
symmetric spacetime, encompassing a vast class of solutions, which when applied to three 
such instances in alternative theories for flavour oscillation and two alternative theories 
for spin oscillation put bounds on the parameters of these theories. Implications are also 
discussed. 


* sumanta@iucaa.in 
1 sumantac.physics@gmail.com 


1 



Contents 


1 Introduction 3 

2 A Brief Introduction to Neutrino Spin Oscillation in Gravitational Field 5 

3 Neutrino Spin Oscillation in a Static Spacetime With Spherical Symmetry 6 

4 Neutrino Flavour Oscillation in Some Classes of Alternative Gravity Theories 9 

4.1 Dilaton Induced Gravity Theory. 11 

4.2 Quadratic Gravity Theory. 15 

4.3 Regular Black holes . 18 

5 Neutrino Helicity Flip in Alternative Gravity Theories 20 

5.1 Einstein-Maxwell-Gauss-Bonnet Gravity. 22 

5.2 f(R) Gravity Theory. 24 

6 Discussion 25 

A Detailed Expressions for Various Quantities 28 

A.l Expressions Related to Spin Oscillation Frequency. 28 

A.2 Expressions Related to Flavour Oscillation Frequency . 30 

B Circular Orbit in Spherically Symmetric Spacetime 31 


2 









1 Introduction 


Neutrino oscillation serves as one of the most cultivating fields in the discipline of elementary 
particle physics. Research along these lines have been boosted after various new experiments 
on solar and reactor neutrinos started yielding tantalizing results [1-8]. Despite of its direct 
connection to the elementary particle physics neutrino oscillation has important contributions 
to cosmology and astrophysics as well. Even though the neutrino oscillation was put to firm 
experimental grounds only recently, the theoretical prediction comes much earlier from the work 
of Pontecorvo [9] and then subsequently it was generalized for propagation within varying density 
medium in [10, 11]. The oscillation we have discussed so far corresponds to neutrino flavour 
oscillation. There is another oscillation between various helicity states within a single flavour, 
known as neutrino spin oscillation. We could also have a mixture of both flavour and spin 
oscillations, which corresponds to neutrino spin-flavour oscillation. 

One of the key research area in neutrino oscillation is the determination of neutrino mass 
and their mixing angles. This focuses on the important fact that even though the mass squared 
difference between various neutrino flavours and their mixing angles are well known, there is 
ambiguity in the determination of absolute values of neutrino masses. This leads to the famous 
paradigm that neutrino masses are either hierarchical or quasi-degenerate in their very nature. 
These results lend themselves to various important questions, which have been key research 
topics in recent years. These include: non-zero value for the mixing angle 0i3, possible aspects 
of CP violation in neutrino oscillation and importantly the non-vanishing 1 — 2 mixing angle. A 
large number of theoretical models exist to explain the above results. These theoretical models 
include, neutrino mass considered to be degenerate by some sea-saw mechanism, assumption of 
large mixing angle for solar as well as for atmospheric neutrino and use of renormalization group 
equation [12-19]. 

Another important place to look for effects of neutrino oscillation is astrophysics. The effects 
being numerous. Since neutrinos are sterile to most of the interactions, after creation they come 
out unimpeded and contain important information about the source. This is the main reason 
for setting detectors like Icecube in order to detect ultra-high energy neutrino from galactic 
centers and active galactic nuclei. There are also other signatures of neutrino oscillation like 
in the Pulsar kick mechanism devised from spin-flavour oscillation of neutrino. The behavior 
of neutrino in high magnetic field existing surrounding a pulsar can lead to resonant oscillation 
and neutrinosphere. These neutrinospheres can explain high proper motion of the pulsar with 
respect to the neighbouring stars [20,21]. 

Neutrino oscillation was first formulated in flat spacetime which was subsequently generalized 
to curved spacetime [22-25] and the formalism can also be used to test the equivalence principle 
[26]. It was well known that for black hole solutions of general relativity the oscillation phase 
along timelike geodesic yields a factor of 2 in comparison to null geodesics. The reason behind 
being the fact that even though the neutrino moves along null geodesic i.e. ds^ = 0, the neutrino 
being massive satisfies the relation, p^ = —m? [27-29]. This discrepancy of factor 2 between 
oscillation phase along null and timelike geodesic persists for a spherically symmetric metric with 
gtt = [32] and has been reproduced here for the most general of such situations. 

Along with neutrino flavour oscillation, the spin oscillation of neutrino is another important 
aspect. In gravitational field the spin of a particle precess and thus there is a non-zero probability 
of neutrino spin flip. There are two possible ways in which the spin oscillation can be studied. The 
first one corresponds to starting with Dirac Hamiltonian in gravitational field and then obtaining 
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mixing terms among various chirality and thus establishing spin oscillation [33,34]. The other 
option is to work with dynamics of spinning particle in gravitational field following [35,36]. The 
spin oscillation probability for circular motion of neutrino in Schwarzschild spacetime has been 
derived in [37] and the generalization to rotating black holes in general relativity has been done 
in [38]. 

However all the above studies of neutrino spin and flavour oscillations have been performed in 
general relativity (except in [32]). Then it is important to understand how these ideas reconcile 
with alternative theories of gravity since there is a general belief that general relativity is only 
a low energy approximation of an underlying fundamental theory [39,40]. In this spirit the 
Einstein-Hilbert action for general relativity gets modified by addition of higher derivative and 
higher order curvature terms. There exist large number of ways in which these higher order 
terms can be introduced in the standard Einstein-Hilbert action. However the criteria that field 
equations should remain second order in the dynamical variable (otherwise some ghost fields 
would appear) uniquely hxes the action to be the Lanczos-Lovelock action [41-46] . Another such 
model explaining the above mentioned problems is obtained by replacing R, the scalar curvature 
in the Einstein-Hilbert action by some arbitrary function of the scalar curvature f{R)- This 
alternative theory for describing gravitational interaction is very interesting in its on spite, for it 
can provide an explanation to various problems in existing general relativity as well as it confronts 
equally well with various tests of general relativity [47-51] . In this context we should also mention 
dilaton couplings in string inspired models, where non-zero electromagnetic fields gets coupled to 
dilaton [52-54] and results in substantial modihcations to Reissner-Nordstrom solution in general 
relativity. Also we can modify the gravity theory by introducing quadratic curvatures of various 
forms. This has been discussed along with possible modifications of the Schwarzschild solution 
in [55,56]. There also exist black hole solutions, which are regular and physical singularity is 
non- existent. These solutions involve modified gravity, properly dressed additional fields, such 
that energy conditions leading to physical singularity are violated [57]. The first such regular 
black hole solution was obtained by Bardeen with magnetic charge and its relation to non-linear 
electrodynamics was shown by Ayon-Beato and Garcia [58] (henceforth referred to as ABG). 
However the solution was non-exact. The exact solution to Einstein’s equation with non linear 
electrodynamics was obtained by Ayon-Beato and Garcia later in [59-61]. We have applied the 
formalism devised by us for the most general static spherically symmetric spacetime to these 
three black hole solutions described earlier and have studied the effect of these alternative theories 
to the neutrino flavour oscillation. 

Eor neutrino spin oscillation we have considered two more alternative gravity theories both 
having interesting theoretical properties. The first one corresponds to vacuum solution in f{R) 
gravity. As we have already mentioned f{R) gravity is a very interesting alternative theory as 
it can explain all the three cosmological phases (for an alternative model see [62,63]), it passes 
through various local tests for general relativity [47-51], moreover it can explain absence of gravi¬ 
ton Kaluza Klein modes in LHC [64] . Motivated by these successes we consider vacuum solution 
in /(i?) gravity which corresponds to the Schwarzschild (Anti) de-Sitter solution. Secondly, we 
consider the second order Lanczos-Lovelock term with additional Maxwell field in addition to 
the Einstein-Hilbert Lagrangian. This is known as the Einstein-Maxwell-Gauss-Bonnet (EMGB) 
gravity. This alternative theory as well posses static spherically symmetric solutions. We have 
considered the spin oscillation probability in these two alternative theories and imposed bounds 
on the parameters. It turns out that these bounds are consistent with previous results. 

In brief, in this work we have generalized previous works on neutrino spin and flavour os- 
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cillation to most general static spherically symmetric spacetime. To our surprise even in this 
general context we can make interesting predictions like, for a neutrino on a circular orbit in 
high energy limit the oscillation frequency vanishes. Hence a high energy neutrino will not have 
any spin flip as long as it is on a circular orbit. We have also derived the oscillation frequency 
for geodesic motion and have used it to compute the probability of helicity flip. On the other 
hand we have also derived oscillation phase and oscillation length for neutrino flavour oscillation 
which has been used to constrain the alternative theories. In this context we have an interesting 
result that oscillation length depends only on the gtt component. It does not depend on any 
other metric components except for conserved energy and angular momentum. 

The paper is organized as follows: we start with a broad introduction to the spin oscillation 
in presence of gravitational held in Sec. (2) and then to spin oscillation frequency in Sec. (3). 
After that we consider general formulation for neutrino havour oscillation and the application 
to alternative theories in Sec. (4). Finally, we provide our analysis for neutrino helicity hip in 
Sec. (5) with concluding remarks in Sec. (6). The detailed calculations are being provided in the 
two appendices, Appendix A, which contain all relevant calculations for both spin and havour 
oscillation and in Appendix B, where we have illustrated the circular orbit related issues for our 
general static, spherically symmetric spacetime. 


2 A Brief Introduction to Neutrino Spin Oscillation in 
Gravitational Field 


Spin angular momentum of a body gets affected by the curvature of spacetime, for example, 
a gyroscope orbiting around a massive body will undergo a precession of its spin. Motion of 
a spinning particle in gravitational held is discussed with great detail in Ref. [35]. Spin of a 
spinning particle can be presented by the use of the spin tensor and it’s momentum such 
that we can introduce the spin vector 

Sp = ( 1 ) 


where, we have introduced the completely antisymmetric tensor density the determinant 

of the metric g = det(p^i,) and the on-mass shell condition, = —mf. For point particles, 

using principle of general covariance it can be shown that the evolution of spin vector and 
four-velocity are parallel transported along its world-line. Also the fact that spin and velocity 
four-vectors are orthogonal to each other remains true even in curved spacetime i.e. = 0. 

However, the properties of the particle are determined by the spin vector in it’s rest frame. 
In order to use this property we should write the geodesic equation in a locally inertial frame. 
This can be done with the help of vierbein vectors e^p '^, where (a) is the Minkowski index. Thus 
we have the following properties of the vierbein vector 




XU _ pU ( 


). x(a) _ p{o.) U . 


( 6 )> 


V{a){b) — Ui' 


( 2 ) 


where we have the following relation: for the inverse vierbein and 'r]{a){b) = 

diag (—1,1,1,1) is the metric tensor in the local Minkowski frame. Thus with the use of vierbein 
vectors the components of spin and four-velocity in a local inertial frame takes the following 
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form, along with the evolution equation 


g{a) ^ 


y(a) 




7 


S(b) 


ds(“ 

dt 




(3) 

(4) 


where, 7 = 1 /° = (dt/dr) and is like the elec¬ 

tromagnetic field tensor and the object, 7(a)(&)(c) = d(a)(d)eJ‘j,^e^^^represents the Ricci 
rotation coefficients. In order to make the spin evolution in the particle’s rest frame, we should 
apply boost within a locally inertial frame. For that we can introduce two four-vectors i?(i) and 
B(,) constructed from such that, and 

Then both “electric” field and “magnetic” field i3(a) has only spatial components denoted 
by E and B such that, = G(o){i) G(i)(j) = ^(i)(i)(/c)^(fe)- These two vectors will govern 

the spin evolution in the rest frame of the particle, which is a linear equation in the spin vector. 
This spin evolution is similar in spirit to the evolution of a charged particle interacting with 
external electromagnetic field. 

Now we can consider three-dimensional spin vector in the rest frame of the particle given by 
Evolution of this three dimensional vector is determined by the following equation: 


f = - (C X G) 
dt 7 

where the vector G is constructed out of two three-vectors E and B as: 


G = 


B 


{E X u} 


(5) 


( 6 ) 


Then the neutrino spin precession is being determined by the vector SI = G/7 and is intimately 
connected to the spin oscillation frequency of the neutrino. 

Hence the main framework to determine the spin precession vector goes as follows: (a) 
Given a spacetime, we need to evaluate the vierbein vectors and the four velocity of the neutrino 
traveling in this spacetime, (b) Then we need to calculate the covariant derivative in order to 
determine the tensor (c) Starting from the tensor calculate the “electric” and 

“magnetic” held and hnally (d) Evaluate the three-vector G in the vierbein frame to get the 
oscillation frequency. The above calculation has been performed in Appendix A.l and we will 
apply these results in the next section to obtain spin oscillation frequency for the most general 
static spherically symmetric spacetime. 


3 Neutrino Spin Oscillation in a Static Spacetime With 
Spherical Symmetry 

The basic equations governing the spin of a particle in a gravitational held has been discussed 
in Sec. (2). In this section, we will apply those equations to determine the spin evolution 
and oscillation of neutrino in a spherically symmetric, static spacetime. This will enable us to 
determine properties connected to neutrino spin as it propagates around a black hole or a massive 
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object. This question was addressed in the context of general relativity, however in this work 
we will try to observe the effect of additional correction terms to the Einstein-Hilbert action 
on the spin oscillation of neutrino and possible shift of oscillation frequency. In this section we 
will outline the basic equations necessary for determination of spin oscillation and then shall 
concentrate on neutrino flavour oscillation in various alternative theories of gravity before again 
returning to neutrino spin oscillation. 

With possible applications in mind the static spherically symmetric metric ansatz has been 
chosen in the form: 

= —f{r)dt^ + (7) 

The horizons in the spacetime are determined by the two conditions: /(r) = 0 and g{r) = 0. 
Except for some simple configuration where f = g, the two conditions in general do not match. 
If we choose our coordinates such that, one of them, namely r, remains constant on the horizon, 
then the normal, oc VqT, has norm = g{r). Thus the condition g(r) = 0 identifies the 
horizon as a null surface. Similarly the condition f(r) = 0, identifies the respective surface as 
a surface of infinite redshift and is known as the ergo surface. Eor spacetime admitting Killing 
vector for time translation symmetry, the norm of the Killing vector = 0 implies that 
/(r) = 0 determines the Killing horizon. Hence both these surfaces have physical significance. 

The oscillation frequency can be obtained in the following steps: First, we need to determine 
the 4-velocity of the neutrino, in terms of the conserved energy and angular momentum and the 
vierbein vectors. Then the covariant derivatives of the vierbein vectors are related directly to the 
“electric” and “magnetic” fields. Thus by evaluating the covariant derivatives we have obtained 
the components of “electric” and “magnetic” fields. Then we can construct the vector G and 
thus the spin precession angular velocity vector SI. As we are in a spacetime which is spherically 
symmetric, we can exploit this symmetry to confine our discussion to the equatorial i.e. 9 = (7r/2) 
plane. There only a single component of spin precession survives. This component has distinct 
expression in different situations. Below we discuss possible situations with respective expressions 
for the spin precession H 2 [for detailed calculation see Appendix A] : 

• We start our analysis with a situation, where the neutrino is assumed to be massive and is 
moving in a circular orbit of radius tq around the black hole. Any orbit around a massive 
object in spherically symmetric spacetime has two constants of motion, namely, the energy 
£ and the angular momentum C. However for a massive particle it is convenient to work 
with energy and angular momentum per unit mass, i.e., f/m and Cjm. In the discussion 
below by energy and angular momentum for a massive particle, will always imply energy 
and angular momentum per unit mass. 

Given the radius rg of the circular orbit we can determine explicitly the energy per unit 
mass E and angular momentum per unit mass L in terms of rg. The expression for four 
velocity, energy and angular momentum so obtained is given explicitly in Appendix B. Then 
we can determine the spin precession frequency explicitly with the following expression: 


1 jW _[[[of_ 

V27(f27+V2/-ro/') 


( 8 ) 


Note that even though we have two functions / and g, the spin precession is mostly 
dominated from the contributions of /, while g appears for a ride. Note that in the limit 
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oi g = f = — {2M/r), we obtain: 



which was derived earlier [37] for circular orbit in the Schwarzschild spacetime and here 
verified from our general analysis. Thus we have derived the expression for spin precession 
frequency for a massive neutrino moving in a circular orbit in the spherically symmetric 
spacetime of Eq. (7). 


• In the literature quite frequently neutrino is taken to move along null trajectory. Thus it is 
important to ask, what happens to the spin precession as the neutrino is taken to move on 
photon circular orbits. This can be obtained directly from Eq. (8) by taking appropriate 
null limit. For circular orbit the null limit is approached when one approaches the photon 
circular orbit, i.e., circular null geodesics located at rg = (2///') (see Eq. (96) in Appendix 
B). Note that since energy and angular momentum are defined as respective quantities per 
unit mass the null limit works out finely. It turns out that by taking appropriate limit i.e., 
To —>■ (2///'), the spin oscillation frequency identically vanishes, leading to. 


n: 


2 (null) 


lim 

■>( 2 ///') 


fW f 

2/ ro 

V 2ro 2/ 

-1 

0 

1 

1 —1 

> 

+ 

1 —1 

_1 


= 0 


( 10 ) 


The same result can also be argued from the expression for Q 2 obtained explicitly for null 
trajectory. Following Appendix B, it turns out that V ,2 scales as 1/E, where E is the energy 
of the neutrino on the null trajectory, which for high energy limit, would tend to the above 
null result. Note that this result is true for arbitrary static and stationary spacetime. For 
a very wide class of solutions this result remains true, as we have not used any specific 
form of the metric keeping the spacetime completely general. 


• In general, if we consider an orbit with non-zero energy and angular momentum then we 
have an expression for spin oscillation frequency in general. It turns out that for zero 
angular momentum the oscillation frequency vanishes. Hence for radial motion there is 
no oscillation in the spin space and thus the spin of the neutrino remains fixed for such 
trajectories. Hence for a general timelike orbit the oscillation frequency takes the following 
form. 


Qgeod _ 's/d 

2 ~ 2Ar2 


1 - 


fEr 


2/(A+ 77) 


( 11 ) 


In the above expression E stands for the energy per unit mass and L accounts for the 
angular momentum per unit mass of the particle moving on the geodesic. Since the mass 
of neutrino is very small, it is highly relativistic. Thus we are really interested in the high 
energy limit of the above expression, which can be obtained by assuming E ^ L and 
then neglecting higher order terms of the ratio L/E. This finally leads to the following 
expression for spin oscillation frequency: 




Lfs/g 

2Er^ 

8 


fr 

2 /. 


( 12 ) 















From the above expression it is evident that for r = 2f/f', vanishes in the high 

energy limit. From Eq. (96) we find that this is precisely the photon circular orbit. Hence 
our previous result can directly be verihed from this line of arguments as well. 

We have derived the spin oscillation frequency of a neutrino along various trajectories. Among 
them we have shown that our general result for spin oscillation frequency in circular orbits 
matches with previous results in the limit of Schwarzschild spacetime. Also we have demonstrated 
that the oscillation frequency vanishes along null circular trajectories. This was shown earlier in 
the context of Schwarzschild spacetime, while in this case we have shown the applicability of the 
result from a more general standpoint. In later sections we shall apply these results to determine 
helicity flip of neutrino and associated parameter space for its detection. 


4 Neutrino Flavour Oscillation in Some Classes of Alter¬ 
native Gravity Theories 


In this section, we will start by first reviewing some basic results about two flavour neutrino 
oscillation in flat spacetime. The eigenstates of the neutrino in flavour basis is denoted by 
|i^flavour)) which is considered as a superposition of eigenstates in mass basis, denoted by 
Then the transformation from flavour basis to mass basis is considered to be performed by 
an unitary transformation U and an associated phase exp(—such that, — 

km ■ Here, Ek^ and km represents the energy and momentum of the fc-th mass eigenstate. 
Then for motion of a neutrino from a point A to another point B the phase in general can be 
interpreted as, 

^km=- r (13) 

J A. 


where, km signifies the fc-th mass eigenstate with ruk being the corresponding mass. This is 
considered as generalization of neutrino phase to curved spacetime if we interpret the four- 
momentum of the neutrino as, = TOfcp^i,(dx'^/(is), where, g^i, is the spacetime metric and s 

is the proper time along the trajectory of the particle. Then it turns out that in the high energy 
limit i.e. Ek S> rrik along with weak held approximation we arrive at the standard oscillation 
phase for two Favour neutrino oscillation [24]. 

Let us start our analysis by hrst considering neutrino Favour oscillation in a general static 
spherically symmetric spacetime. Then after deriving the expressions for oscillation probability 
and oscillation length in the general case, we will apply these results for some speciFc situations 
in alternative theories, which in turn will constrain various parameters of these models. For 
these purposes we will follow Ref. [32] . The detailed expressions are being provided in Appendix 
A, we will summarize the important results and shall provide the physical insights behind these 
results. 

For the general metric ansatz presented in Eq. (57) we have two conserved quantities the 
energy and angular momentum. Also the motion is taken to be conFned in the equatorial plane. 
Thus the phase along null geodesic turns out to be: 





rrik // 

2EkmVv\l g_ 


(14) 
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where, the potential V turns out to have the expression: V = 1 — (fLl/r^El) and the index 


k merely specifies that we are considering kth massive neutrino. The standard result for phase 
can be obtained from the high energy limit, with V ^ 1 and f = g = (1 — 2M/r). This leads to 
the following expression for the null phase: 



(15) 


which exactly matches with the respective result for Schwarzschild spacetime [24], As the neu¬ 
trino travels with speed very close to that of light it is considered to have traveled along a null 
line. However it would be more appropriate to calculate the phase along a timelike trajectory. 
These two phases differ necessarily as we have taken neutrino to be massive but moving along 
null geodesic. Then for timelike geodesic the phase turns out to be: 



(16) 


Thus we found that phase along null and timelike trajectory differ by a factor of 2 even for 
the most general spherically symmetric spacetime discussed here. This factor exists in the 
flat [28], Schwarzschild [25,27], Kerr-Newmann [30] and in various alternative gravity theories 
with solutions having gu = g'"^ [32]. If we neglect the coherence effects, i.e., assume that 
neutrinos of different flavours are created at the same spacetime point and are also detected 
at the same spacetime point, then this phase difference is irrelevant for interference. In this 
particular situation as well we have approximated the trajectory of ultra-relativistic neutrinos 
as null geodesics which starts and ends at the same spacetime point and thus this factor would 
not show up in interference. However in principle one could have chosen more complicated 
trajectories as well (see for example [65]). 

However as argued in [66] this factor appears in flat spacetime due to wrong use of the 
group velocity in phase. In curved spacetime as well its origin is embedded in the following 
reason—we are taking the neutrino to be massive, i.e., PaP“ = but at the same time we 

are approximating it to be moving on a null geodesic since it is ultra relativistic. Due to this 
erroneous use of two separate formulations in a single problem we are getting this factor of 2. 

Having obtained the oscillation phase it is important to consider some non-trivial effect 
that these alternative theories have on top of the standard Schwarzschild coordinates. The best 
physical observable for this purpose as pointed out in [32] is the oscillation length. The oscillation 
length can be obtained as [67] 



(17) 


Then the differential form of the oscillation phase can be written in terms of differential of proper 
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length from Eq. (14) as 


d<i) 


null 

km 


rukdr // 

2Ek^Vv\l 9 


mkVf 


dl 


2 

mi 


Vfdl 


(18) 


The remarkable fact about this equation is that the differential phase does not depend on the 
Qrr component. Thus only the gtt component appears in the oscillation phase and thus affect our 
solutions. There is another assumption that has gone into the above equation, which corresponds 
to taking the eigenstates of mass and energy identical. This has been scrutinized critically by 
several authors and in general it is assumed that pg ^ represents common energy of mass 
eigenstates. Also there is one critical point that we should stress, as the flat spacetime has both 
time and space invariance, both pg and p^ carry equal momenta. However in curved spacetime pg 
is a conserved quantity due to existence of timelike Killing vector held, which obviously is not true 
for pr- Hence the equal momentum assumption does not work in curved spacetime [25,29,31]. 

Hence the phase shift determining the oscillation turns out to be, d^k — d^j oc dl. 

Then we can dehne the oscillation length Lose for neutrino in curved spacetime as. 


^grav 



47rpg 1 

^/fF) 


(19) 


which for a hat spacetime reduces to, Log* = (^Trpg/Am^^). Hence hnally the fractional change 
in oscillation length in presence of gravity turns out to be 


A/i 




( 20 ) 


while the other object of interest turns out to be the change in oscillation length in alternative 
theories in comparison to the corresponding solution in General Relativity. This one has the 
following expression 


Ah 


1 

\//alt(r) 


1 

\//gr(?’) 


( 21 ) 


In the above expression /ait(?') is the gtt component in the alternative theory and /gr(?') is the 
corresponding gtt element in spacetime described by General Relativity. Next we will try to 
compute these changes to put bounds on parameters in the theory, which we will compare with 
the respective parameters using spin oscillation as well. 


4.1 Dilaton Induced Gravity Theory 

In General Relativity, Schwarzschild solution represents the vacuum solution outside a spherically 
symmetric massive object. If the mass of the central object becomes comparable to Planck scale 
then also the same solution describes quite well the structure of the spacetime, except for regions 
near the massive object (if it forms a black hole, then near the singularity). However when we 
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bring coupling with Maxwell field then the corresponding classical solution, which is the Reissner- 
Nordstrom solution gets modified significantly since every solution with non-zero couples 
with dilaton. Hence the Maxwell Coupled solution gets drastically modified due to presence of 
the dilaton field and thus the low energy four dimensional effective action takes the following 
form as obtained from string theory [52] : 


A = 




-R + -f 2 (V$)^ 


( 22 ) 


In the above expression, — dyA^, is the Maxwell field tensor, belonging to the 17(1) 

subgroup of x Eg,. All other gauge fields except $ have been set to zero, then for a purely 
magnetic Maxwell field we have, the only non-zero component to be, E^^ = Q sin 0, leading to, 
E"^ = 2Q^IR'^. It turns out that only three components of Ricci tensor are non-zero, with one of 
them satisfying the relation, i?oo = —V^d>. Then the spherically symmetric solution turns out 
to be [52]: 


ds^ 



dt^ F 



+ r 





(23) 


where, = d9^ +sm^ OdcjA is the surface element on (d, (j)) plane, $o represents the asymptotic 
value of dilaton and Q, the black hole charge. This metric resembles the Schwarzschild solution, 
however with another horizon being located at r = /M. Both the Reissner-Nordstrom 

and dilaton solution describes a black hole for small Q/M ratio, while in other situations (with 
Q/M < 1) both describe naked singularity. Also from the string theory viewpoint, the strings 
do not couple to the metric g^y but rather to e^^g^y. Then in string a model the effective 
Lagrangian leads to the following action [52]: 


A = 


d^xy/^e-^'^ 


-R-A (V$)^ -f F^yE^^'' 


(24) 


From this modified action we obtain the following solution for the metric in spherically symmetric 
coordinate as [52]: 


ds^ 


1 - (2Me‘*'« /p) 2 V 

1 - ^ [1 - {2Me‘^o/p)] [i _ /Mp)] 


(25) 


Note that this line element is exactly in the form of the metric ansatz we started given in Eq. (7). 
The most important parameter in this theory is the dilaton charge, which can be obtained by 
integrating over a two sphere at spatial infinity as [52,56], 



^ 2M 


(26) 


Having obtained the line element we will now compute the oscillation probability and thus the 
oscillation length. The most important object that appears in all the expressions correspond to 
the potential V{r), which in this particular case turns out to have the expression: 


V{r) = 1 - 


ll^ l-(2Me‘^Vp) 
AEl^ 1- (Q2e3<[-o/Mp) 


(27) 
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Figure 1: (color online) In the two figures we have plotted change in neutrino oscillation length 
induced by dilaton coupled Maxwell Field. The left figure shows variation of Ali difference of 
neutrino oscillation length from the dilaton coupled theory to that of flat spacetime with radial 
distance from the source. Thus as we move outwards the difference goes to zero, since at large 
distance from the source the spacetime becomes flat. While the second figure depicts variation of 
AI 2 the difference from Schwarzschild geometry. Also at large distance this difference approaches 
zero, showing at large scale the solution resembling Schwarzschild spacetime. 


Hence the phase along null geodesic as well as that along timelike geodesic can be computed in a 
straightforward manner by just substituting the above expression for potential and corresponding 
expressions for metric elements in Eq. (14) and Eq. (16) respectively. Then the two oscillation 
length, first one corresponding to difference from flat spacetime have the following expression: 


Ah 


l_(Q2g3$o/Mp) 
1 - i2Me^o/p) 


(28) 


This expression has very little dependence on the charge Q since for normal astrophysical systems 
Q ^ M. The other comparative oscillation length corresponds to the difference from the general 
relativistic counterpart i.e. the Reissner-Nordstrom solution. This comparative oscillation length 
turns out to have the following expression: 


1- (Q2g3<I.o/Mp) 
l-(2Me‘^o/p) 




Finally, we can also compute the oscillation probability for an electron type neutrino to remain 
an electron type neutrino for various lengths. Then we can compare it with the solar neutrino 
result in order to constrain various parameters of this theory. 

The variation of the two comparative oscillation lengths with radial distance have been pre¬ 
sented in Fig. 1. Also the probability for an electron type neutrino to remain another electron 
type neutrino for a definite energy window has been presented in Fig. 2. The fact that we 
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Figure 2: (color online) In this figure we have depicted the probability of an electron type neutrino 
to remain an electron type neutrino with its energy in MeV for a length of 180 km. Different 
curves describe the transition probability corresponding to different values of , which depends 
on the asymptotic value of the dilaton. Green curve represents oscillation probability in absence 
of dilaton field, while the other two curves depict oscillation probability with dilaton field being 
present. 
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Table 1: Results from real time experiments regarding solar neutrino flux have 
been shown. The errors presented are statistical errors. Bounds on the dilaton 
charge D from each of these experiments have been estimated. 


Experiment 

Reaction 

V flux 

Bound on dilaton 
charge D 

Kamiokande [1] 

ve 

2.80 ±0.19 

< 

2.12 

X 

10- 


Super-K I [2] 

ve 

2.38 ±0.02 

< 

2.05 

X 

10- 

-9 

Super-K II [3] 

ve 

2.41 ±0.05 

< 

2.07 

X 

10- 

-9 

Super-K III [4] 

ve 

2.32 ±0.04 

< 

1.94 

X 

10 - 

-9 

SNO Phase I [5] 

cc 

1 7f)+0-06 

-L- <0_Q Q5 

< 

1.74 

X 

10 - 

-9 

(pure D 2 O) 

ve 

2 QQ+0.24 

z.ot;_Q 23 

< 

2.05 

X 

10 - 

-9 


NO 

5.09tn3 

< 

3.76 

X 

10 - 

-9 

SNO Phase II [6] 

CC 

1.68 ±0.06 

< 

1.69 

X 

10 - 

-9 

(NaCl in D 2 O) 

ve 

2.35 ±0.22 

< 

1.98 

X 

10 - 

-9 


NC 

4.94 ±0.21 

< 

3.53 

X 

10 - 

-9 

SNO Phase III [7] 

CC 

-0.04 

< 

1.68 

X 

10 - 

-9 

{^He counters) 

ve 

1 '7'7+0.24 
^-0.21 

< 

1.74 

X 

10 - 

-9 


NC 

r r^+0.33 
0-04_0.31 

< 

3.92 

X 

10 - 

-9 

Borexino [8] 

ve 

2.4 ±0.4 

< 

2.07 

X 

10 - 

-9 


are considering probability of an electron type neutrino to remain an electron type neutrino is 
illustrated in Fig. 2 as the maximum probability reaches the value unity. 

However we can do more, from the neutrino oscillation probability, in this gravity model 
we can constrain the parameter D, defined through Eq. (26) using the data for solar neutrino 
oscillation. This can be performed along the following lines, since the geometry and hence 
the oscillation probabilities are different, the theoretical flux would differ from the Schwarzschild 
value and this difference depends upon the dilaton charge. To be compatible with the experiment 
the excess flux must be within the statistical errors present in observed fluxes. This in turn will 
provide stringent constraints on the dilaton charge D. The constraints on the dilaton charge for 
various solar neutrino experiments have been presented in Table 1. 

4.2 Quadratic Gravity Theory 

In high energy physics general relativity is often treated as the low energy approximation of 
some underlying fundamental theory. This idea is reconciled by adding higher curvature terms 
to general relativity and observing their implications in high energy phenomenon through astro- 
physical experiments. For example, modification of accretion disk structure in presence of higher 
curvature term has been discussed in Ref. [71]. In this section we consider a class of alternative 
theories in four dimensions, which is obtained by modifying the Einstein-Hilbert action through 
introduction of various quadratic and algebraic curvature scalars with proper couplings, such 
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that the action gets modified to [55,56,68]: 


A = d^Xy/^ 


R 


167rG' 


aifi{v)R^ + a2f2{v)Rf,uR^'' + azf2,{v)R^yapR^''°'^ 


+ aiU{v)R%^pR^^''‘^^ - ^ {Vav^^^v + 2V{v)} + 


(30) 


where g represents the determinant of the metric g^i,. Among other quantities we have the Ricci 
scalar R, Ricci tensor Riemann tensor R^vap along with its dual constructed from 

the metric gap- In the above action Amatter represents the matter Lagrangian, v is an arbitrary 
scalar field with {ai,(3) representing coupling constants. These theories have their motivation 
in the low energy expansion of string theory [69,70]. Then by varying the above action field 
equation with respect to g^i, can be obtained, to which spherically symmetric solution can be 
arrived at through the following metric ansatz 


- /o [1 + <^do{r)] dt^ + /g ^ [1 + efco(r)] dr'^ + (31) 


where e comes from expansion of the scalar field around its solution and /o = 1 — (2Mo/r). Here 
Mq represents the bare mass and dfl^ is the two dimensional surface element. Also in the above 
equation /ig and fcg represents small deformation around the Schwrazschild value. Then we can 
use solution of scalar field equation in order to obtain a solution to the modified field equation at 
linear order in e. Further the requirement that metric should be asymptotically flat and regular 
at r = 2Mq can lead to unique solution. Then defining a dimensionless coupling parameter, ( = 
{lenGal/pMi), we can introduce a physical mass parameter given by M = Mq [1 + (49/80)^]. 
Then the modified line element takes the following form [55,56]: 


ds^ 


-f{r) 


1 + 


3/(0 



dR 


fir) 




dr^ + r^dn^ 


(32) 


where we have /(r) = 1 — (2M/r), and the two other unknown functions h{r) and fc(r) have the 
following expressions: 


26M 66 96 80M‘^ 


M0-1 + —+ y^ + y^- 


,, , ^ M 52 M2 

k(r) = l + - + - — 


2M3 16M4 


5A 


368 M^ 
3 r® 


(33) 

(34) 


The important point to be stressed is that the metric elements so obtained do not depend on 
the bare mass Mq but on the physical mass M. Thus having obtained the solutions we will now 
consider the oscillation length and oscillation phase of neutrino in this spacetime. 

The line element is exactly in the form of the metric ansatz we started with as given in 
Eq. (7). Having arrived at the line element we will now compute the oscillation probability and 
thus the oscillation length. The important expression that appears in all the expressions for 
oscillation phase is the potential V (r), which in this case leads to the following expression: 


V{r) = 1 - 


/2 


r^El 


-fir) 


1 + 


3fir) 



(35) 
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Figure 3: (color online) In this figure the variation of comparative oscillation length A^i and 
AI 2 with radial coordinate have been presented for various choice of parameters. The left figure 
shows variation of A^i, the difference of oscillation length in quadratic gravity and fiat spacetime. 
On the other hand, the second figure depicts variation of A/ 2 , difference from Schwarzschild 
geometry. At large distances this comparative difference tends to zero resembling Schwarzschild 
behavior. 



Figure 4: (color online) In this figure we have depicted the probability of an electron type 
neutrino to remain an electron type neutrino with its energy in MeV for a length of 180 km. 
Different curves describe the probability of this event corresponding to different values of 
Green curve represents oscillation probability in absence of quadratic corrections, while the other 
two curves depict oscillation probability with quadratic corrections being present. 
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where, f{r) = l — {2M/r) and h{r) is given by Eq. (33). Then the phase along null geodesic as well 
as that along timelike geodesic can be computed directly by substituting the above expression for 
potential and corresponding expressions for metric elements in Eq. (14) and Eq. (16) respectively. 
Then the two oscillation length, first one related to the difference from flat spacetime have the 
following expression: 


Ah 


fir) l + %) 


- 1 


(36) 


while the other one corresponding to the difference from the general relativistic counterpart i.e. 
the Reissner-Nordstrom solution turns out to be 


Ah 


fir) 


1 + 


3/(r 


iff Hr) 



(37) 


In both the expressions, /(r) = l — {2M/r) and h{r) is given by Eq. (33). Finally, we can compute 
the oscillation probability of electron type neutrino converting to electron type neutrino which we 
can compare with the solar neutrino result in order to constrain various parameters of this theory. 
Even though the parameter ^ appears explicitly in the metric elements, it is not a fundamental 
parameter. The fundamental parameter is which is related to C via C = {I^ttGo^Mq). 
Figure 3 depicts the difference lengths Ah and Ah for various choice of parameters, while figure 
4 depicts oscillation probability for various choices of af. In Table 2 we present constraints on 
from solar neutrino experiments using identical techniques as explained in previous section. 


4.3 Regular Black holes 


Existence of singularity appears to be an inherent property of most of the solutions to gravita¬ 
tional field equation in general relativity. This problem is generally avoided by considering cosmic 
censorship conjecture according to which singularities are always dressed by event horizons. This 
shows that any pathological behavior at the singularity has no influence on the exterior region. 
To circumvent these difficulties a set of regular black hole solutions were proposed, known as 
“Bardeen black holes” [72]. However none of these models is an exact solution to Einstein’s 
equation with some known physical source associated. It was first suggested in [59-61] that 
even within the context of general relativity it is possible to construct singularity free solutions. 
However this can be achieved only at the price of introducing non-linear sources. Thus by in¬ 
troducing non-linear electrodynamics to Einstein gravity it is possible to obtain singularity free 
solutions. The solution so obtained has the line element with /(r) = g{r), where the function 
/(r) has the following expression: 


1. Bardeen spacetime 


2. ABG spacetime 


fir) = 1 - 


2mr^ 

(^2 +^ 2 ) 3/2 


2mr^ 

(^2 -(-^2)3/2 (g2_j_^2^2 


(38) 


(39) 
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Table 2: Results from real time experiments regarding solar neutrino flux have 
been shown. The errors presented are statistical errors. Bounds on from each 
of these experiments have been estimated. 


Experiment 

Reaction 

V flux 

Bound 

on 


Kamiokande [1] 

ue 

2.80 ±0.19 

< 

0.95 

X 

10 - 

is 

Super-K I [2] 

ue 

2.38 ±0.02 

< 

0.42 

X 

10“ 

18 

Super-K II [3] 

ve 

2.41 ±0.05 

< 

0.45 

X 

10 - 

18 

Super-K III [4] 

ve 

2.32 ±0.04 

< 

0.36 

X 

10 - 

18 

SNO Phase I [5] 

GO 

-L- Q5 

< 

9.16 

X 

10 - 

19 

(pure D 2 O) 

ue 

^•oy_o.23 

< 

0.43 

X 

10 - 

18 


NC 

o.uy_o 43 

< 

2.32 

X 

10 - 

18 

SNO Phase II [6] 

GG 

1.68 ± 0.06 

< 

8.23 

X 

10 - 

19 

(NaCl in D 2 O) 

ve 

2.35 ±0.22 

< 

0.39 

X 

10 - 

18 


NC 

4.94 ±0.21 

< 

2.16 

X 

10 - 

18 

SNO Phase III [7] 

cc 

-0.04 

< 

8.12 

X 

10 - 

19 

{^He counters) 

ue 

1 7'7+0.24 
'- 0.21 

< 

9.19 

X 

10 - 

19 


NC 

c r4+0.33 
^•^^-0.31 

< 

2.58 

X 

10 - 

18 

Borexino [8] 

ve 

2.4 ± 0.4 

< 

0.44 

X 

10 - 

18 


In the above expressions m denotes the standard gravitational mass and q stands for charge 
parameter measured in units of mass m. Since both the line elements are exactly in the form 
presented in Eq. (7) we can carry forward our analysis presented in earlier sections. The impor¬ 
tant expression appearing in all the related expressions is the potential term V(r) having the 
expression: 

1. Bardeen spacetime 




1 - 


2mr^ 


(q2 +^2)3/2^ 


2. ABG spacetime 


V{r) = 1- 


r^El 


1 - 


2mr^ 


+ 


9 2 


(g2 _|_ ^2)3/2 


(40) 


(41) 


Then the above expressions when substituted in Eqs. (14) and (16) the respective expressions 
for oscillation phase along null and timelike geodesics can be obtained. Also the two oscillation 
lengths, one related to the oscillation length difference between flat spacetime and the curved 
spacetime has the following expression: 
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1. Bardeen spacetime 


Ah 


1 



2mr^ 

(g2+r2)3/2 


- 1 


(42) 


2. ABG spacetime 


Ah 


1 



2mr^ 

(q2_,_^2)3/2 


W+^ 


- 1 


(43) 


The second one represents the difference between the oscillation length in the spacetime governed 
by alternative gravity theories and corresponding general relativity solution. This has the 
following expression: 

1. Bardeen spacetime 


Ah 


1 



2mr^ 

(g2+^2)3/2 



(44) 


2. ABG spacetime 


Ah 


1 



2mr^ 


(q,2_|_^2)3/2 + (q,2_|_^2) 





(45) 


Both the comparative behaviors regrading oscillation lengths have been illustrated. Figure 
5 depicts the comparative oscillation lengths Ah and Ah for Bardeen spacetime and figure 6 
shows identical diagrams but for ABG spacetime. Having obtained all the oscillation lengths 
and their comparative behavior, we can now compute the oscillation probability of electron type 
neutrino converting to electron type neutrino. It turns out that the corrections to the oscillation 
length and hence the departure of oscillation probability in these regular black hole solutions 
from that in the Schwarzschild spacetime is quite small. They cannot be used accurately to 
place tight constraints on the parameter q as we have done in the previous two situations. This 
originates from the fact that all the corrections are sub-leading and falls off much faster, leaving 
very little correction in comparison to the previous two cases. 

Having discussed three alternative gravity theories and the neutrino flavour oscillation thereof, 
we will now concentrate on neutrino helicity flip, i.e., spin oscillation of neutrino in alternative 
theories. 


5 Neutrino Helicity Flip in Alternative Gravity Theories 

In Sec. (3) we have derived the spin oscillation frequency for neutrino, both along circular and 
non-circular geodesics. We have observed that for circular motion in the high energy limit, the 
oscillation frequency vanishes. This is a very important result, it suggests that a neutrino moving 
in a circular orbit does not undergo any change of its spin and hence helicity (since neutrino is 
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Figure 5: (color online) In this figure the variation of comparative oscillation length A^i and 
AI 2 with radial coordinate has been presented for various choices of parameters. The left figure 
shows variation of Ali, the difference of oscillation length in Bardeen spacetime and flat space- 
time. On the other hand the second figure depicts variation of A/ 2 , difference from Reissner- 
Nordstrom geometry. At large distance this comparative difference tends to zero resembling 
Reissner-Nordstrom behavior. 




Figure 6: (color online) In this figure the variation of comparative oscillation length A/i and A /2 
with radial coordinate has been presented for various choices of parameters in ABG spacetime. 
The left figure shows variation of A/i and the right figure depicts variation of A/ 2 . See text for 
more discussions. 
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always a highly relativistic particle). More importantly this result is true not only in Einstein 
gravity but in all other alternative gravity theories, as the result essentially depends only on 
static and spherical symmetry. This result was first derived for Schwarzschild spacetime in [37] 
and in this work is shown to transcend general relativity. 

To start with we assume that the neutrino initially is left-handed, implying its spin vector 
being anti-parallel to the velocity of the particle. From Eqs. (3) and (4) from Sec. (2), along 
with Eqs. (8) to (12) in Sec. (3) it is clear that neutrino spin rotates around the second axis. 
This enables us to construct the effective Hamiltonian for neutrino spin oscillations in the static 
spherically symmetric spacetime which takes the following form: 


HeS 


0 —ifl2 
ifl2 0 


(46) 


Thus the effective Hamiltonian depends solely on the oscillation frequency O 2 , which vanishes 
for neutrino moving on a circular orbit (we only consider the high energy limit in the remaining 
discussion). Then from Eq. (46) the neutrino spin oscillation probability after traversing a 
distance r in time t = r/c turns out to be: 

P{t) = sin^ (H2t) (47) 


Then we can read off the oscillation frequency H 2 from Eq. (12) and then substitute in Eq. (47), 
which ultimately leads to the following expression for spin oscillation probability of neutrino 
for the most general static spherically symmetric spacetime in geodesic with energy per particle 
mass E and corresponding angular momentum L in the high energy limit as: 


P{L,E,t)\geod = sin^ 


LfVg \ fr 
2Er2 [ 2/ 


(48) 


The most striking feature of the above expression is that above all other dependencies the 
oscillation probability depends on the angular momentum L. This clearly depicts the fact that 
if we consider only radial motion, the spin oscillation probability identically vanishes. Thus for 
any static, spherically symmetric spacetime neutrinos traveling along radial direction would not 
suffer helicity flip. If it was left handed originally it will remain left handed for radial motion. 

Note that in all the three cases considered above since /(r) and g{r) are asymptotically flat 
and there is an additional 1/r^ factor in the probability. Thus at large distance (for example on 
earth) the probability for helicity flip will be vanishingly small. Thus in this case we consider 
two more alternative gravity theories namely, (a) The Einstein-Maxwell-Gauss-Bonnet gravity 
and (b) f{R) gravity theory both of them are asymptotically de-Sitter or Anti de-Sitter. Thus 
they can produce significant probability for neutrino helicity flip and we can obtain some bound 
on their parameters to ensure that the helicity flip remains within experimental bounds. We can 
also compare these bounds with previously obtained results in [32,56]. 


5.1 Einstein-Maxwell-Gauss-Bonnet Gravity 

Spacetime having more than four dimensions is an interesting concept. For it can solve for 
various fundamental problems in theoretical physics from cosmological constant to hierarchy 
problem [74-77]. In these extra dimensional scenarios, the spacetime we live in is assumed to be 
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a four dimensional brane, embedded in a higher dimensional bulk. Ordinary matter fields are 
confined in the brane, while gravity can propagate in the bulk as well. Also as we have argued 
earlier the Einstein-Hilbert action is supposed to be a low energy realization of the fundamental 
theory. To keep out the ghost terms it is instructive to modify the gravity action by including the 
second order Lovelock term which is known in the literature as the Gauss Bonnet term. After 
inclusion of the Gauss-Bonnet (GB) term and a Maxwell field along with the Einstein-Hilbert 
term the modified action looks like 


5 


dx^y/^ 


R cx 


+ A") + 


(49) 


In the above expression R, R^^, and R^^a/3 ^re Ricci scalar, Ricci tensor and Riemann tensor 
respectively, is the electromagnetic field tensor and a being the GB coupling coefficient with 
dimension of length squared. The field equations for gravity can be obtained by varying the above 
action with respect to the metric gap and variation of the electromagnetic field tensor would 
lead to electromagnetic field equations respectively. Through this variation we obtain [78-81] 


R 




g^^R - a ^-g^^ [Rap-fsR‘^^'^' - + R^) - ^RR^. + 4i?^„i?° 


AR’^PR^aup - 2R ^P^R,ap^ 


= T, 




= 0 


(50) 


where is the usual stress tensor for electromagnetic field. The important thing to notice is 
that the field equation only contains second order derivatives of the metric, no higher derivatives 
are present. This is expected since Gauss-Bonnet gravity is a subclass of Lovelock gravity, which 
does not contain higher derivative terms of the Riemann tensor. 

We can obtain static spherically symmetric solutions to these field equations having the form 
of Eq. (7). It turns out that these solutions are asymptotically de-Sitter or Anti de-Sitter [78]. 
Also we should mention that the solution obtained from the above field equations will be in higher 
dimensions, i.e., the line element would correspond to: ds^ = —f{r)dt^ + f~^{r)dr‘^ -I- r‘^dVl\^ 
where dHg = d0^-|-sin^ 0i(d0|-l-sin^ 92d9^). However we have been emphasizing that this solution 
should be interpreted from a brane world point of view, with the visible brane characterized by 
03 = constant hypersurface. With this choice the solution reduces to four dimensions with 
(t,r, 01,02) as the set of coordinates. Such spherically symmetric solutions were obtained in 
Ref. [78] and has the particular form with reference to Eq. (7) as 


f(r) = g{r) = K+ — 


1 ±\ 1 + 


8a (to -I- 2a j K j) 8ag^ 


3r® 


(51) 


where K determines the scalar curvature of the spacetime which can take values 0, ±1. However 
in this work we attribute to K the value 1. Then form solar system tests and neutrino oscillation 
experiments [32,56] we can infer stringent bounds on a~^. As we have mentioned that we will 
work at a large distance from the source thus for our study we can make a power series expansion 
of the terms inside the square root in inverse powers of r and arrive at the following result 


/(r) = g{r) = 1 + — + 


m + 2a 


q 

3M 


(52) 
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Figure 7: (color online) In this figure the variation of neutrino oscillation probability with neu¬ 
trino energy has been presented for various choices of angular momentum L and the GB param¬ 
eter a. See text for more discussions. 


Then in the large r limit we have the following expression r/'/2/ ^ 1 — {2ajr^). Using the large 
r limit of both /(r) and g{r) we obtain the neutrino helicity flip probability to be, 

P(Z/, U, a)Igeod = sin ^2E\/2^ 

Note that the helicity flip depends on the angular momentum of the particle L, the energy E 
and the GB coupling parameter a. In figure 7 we have plotted the helicity flip probability with 
the energy of the neutrino. It turns out that for all values of a and L the probability of helicity 
flip is larger for low energy neutrino. While for high energy neutrino the helicity flip probability 
is substantially small. Also as the GB parameter a increases (or equivalently as a~'^ decreases) 
the probability of neutrino helicity flip decreases signihcantly. The low energy solar neutrinos 
have energy of about 0.4 MeV. Then for L ^ 100 and a ~ 10^° we get P ~ 0.02. This is quite 
consistent with observations as well [82], where from the Kamiokande-II data the helicity flip 
probability was obtained as 0.07. Thus the helicity flip probability of neutrino leads to a 
bound on the GB parameter which turns out to be a 9.34 timeslO^. 

5.2 f(R) Gravity Theory 

There exist another way of modihcation of the Einstein-Hilbert action, which is obtained by 
introducing a term f{R) in the general relativity lagrangian, where / is taken to be some 
arbitrary function of the scalar curvature R. In order to get the gravitational field equations we 
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will use the standard method, i.e., we will vary the metric leading to [49,83,84] 

\g^.uf{R) - - g^.af'{R) + V^V./'(i?) = -4^T“ (54) 

The vacuum solution, i.e., solution with T]))]®'**®'' = 0 corresponds to a de-Sitter Schwarzschild 
or Anti de-Sitter Schwarzschild solution for which the Ricci scalar is covariantly constant. This 
constant Ricci scalar in turn corresponds to R^y oc g^j^y. Since □/'(i?) = 0 in the scenario we 
are considering the field equation given in Eq. (54) reduces to the following algebraic equation, 
0 = 2/(i?) — Rf'{R). This immediately suggests that the model f{R) oc satisfy the above 
equation [49]. Hence the de-Sitter or Anti de-Sitter Schwarzschild solution is an exact vacuum 
solution the the f{R) gravity theory with the line element being given by (see Eq. (7)) 


-1 




dr^ - 


(55) 


which immediately shows that /(r) = g{r). Here the minus(plus) sign corresponds to the 
spacetime being (Anti) de-Sitter spacetime, M corresponds to the mass of the black hole and 
Lf \s the length parameter of the (Anti) de-Sitter spacetime. This length can be related to 
the scalar curvature as i? = ±w (where also the plus sign corresponds to de-Sitter spacetime 

and minus sign corresponds to Anti de-Sitter spacetime). Then in the large r limit we have 
rf'/2f ^ 1 ± (L^/r^). Using the large r limit of both /(r) and g{r) we obtain the neutrino 
helicity flip probability to be, 

P(L, A, L/)|geod = sin (^2ELf^ 

The expression for helicity flip probability clearly shows that it depends on the angular mo¬ 
mentum of the particle L, the energy E and the length parameter Lf in f{R) gravity. In figure 
8 we have plotted the helicity flip probability with the neutrino neutrino energy for different 
choices of other two parameters. As in the GB scenario in this case as well it turns out that for 
all values oi Lf and L the probability of helicity flip is larger for low energy neutrinos. While 
for high energy neutrinos the helicity flip probability is substantially smaller. Also as the length 
parameter Lf increases the probability of neutrino helicity flip decreases significantly. The low 
energy solar neutrinos have energy of about 0.4 MeV. Then for L ~ 100 and Lf ^ 10® we get 
P ^ 0.02. This is quite consistent with observations as well [82], where from the Kamiokande-H 
data the helicity flip probability was obtained as ~< 0.07. Thus the helicity flip probability 
of neutrino leads to a bound on the length parameter in f{R) gravity which turns out to be 
Lf ~> 9.87 X 10"^ which is consisten with the bound presented in [32]. 


6 Discussion 

Neutrino oscillation has been extensively studied in flat spacetime and in general relativistic 
solutions. Given the recent boost in the search for alternative gravity theories it is legitimate 
to ask the status of neutrino oscillation in these alternative gravity theories. A first step along 
this direction was taken in [32]. In this work we have generalized their setup by considering the 
most general static spherically symmetric spacetime. Our strategy in this work is as follows: 
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Figure 8: (color online) In this figure the variation of neutrino oscillation probability with neu¬ 
trino energy has been presented for various choices of angular momentum L and the length 
parameter Lf in f{R) gravity model. See text for more discussions. 
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• First we have derived general results for an arbitrary static spherically symmetric spacetime 
in the context of both spin and flavour oscillation. For neutrino spin oscillation we have 
shown that the oscillation frequency is dependent on the angular momentum, energy and 
both gtt and grr components of the metric. It turns that for radial motion the spin 
oscillation frequency (and hence the probability) identically vanishes. Also for circular 
orbits in the high energy limit the oscillation frequency vanishes. Both these results hold 
for any static spherically symmetric spacetime. 

• For neutrino flavour oscillation we have derived both the oscillation phase and the os¬ 
cillation length for our general setup. It turns out that the oscillation length depends 
only on the gtt component. Also starting from the oscillation length we have two more 
quantities, namely departure of the oscillation length from both the flat spacetime and 
the corresponding general relativity solution provides an ideal test to obtain departure of 
alternative theories from their general relativity counterpart. 

• We have discussed three alternative theories to illustrate the results obtained in the con¬ 
text of our general static spherically symmetric spacetime. These three theories include 
dilaton coupled Maxwell held, Einstein-Hilbert action modihed with all possible quadratic 
correction terms and hnally regular black hole solutions. By calculating the oscillation 
probability in all these theories and comparing them with the solar neutrino results we 
can put bounds on the parameters in these models. It turns out that for dilaton coupled 
Maxwell held the dilaton charge has the following bound D <~ 1.69 x 10“® from the 
SNO Phase-II data. For the quadratic gravity theory the bound on Og ^ turns out to be 
Og ^ < 0.43 X 10“^® obtained from SNO phase-I data. While for regular black holes no 
such bound was obtained as for them the corrections are within experimental errors. 

• In the case of neutrino spin oscillation, we have considered two separate gravity theory. 
The f{R) gravity theory and the Gauss-Bonnet gravity theory. In both the theories the 
helicity hip probability depends on neutrino energy and it turns out to be signihcant 
when the neutrino energy is small, while the probability is quite small for high energy 
neutrinos. This feature was observed in both the theories. Moreover using the helicity hip 
probability we obtain the following bounds a ~> 9.34 x 10® for Gauss Bonnet parameter 
and Lf ~> 9.87 x 10^ for length parameter in vacuum f{R) gravity. 

Thus, starting from a general static spherically symmetric metric ansatz we have derived both 
neutrino havour oscillation probability and neutrino spin oscillation probability. After obtaining 
these general results we have applied them to various static spherically symmetric solutions in al¬ 
ternative gravity theories. These in turn when compared with experiments produce experimental 
bounds on various parameters in these alternative theories. 
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A Detailed Expressions for Various Quantities 

In this and subsequent sections we provide a detailed analysis of the expressions. We have not 
included these results in the main text in order to maintain the flow of ideas in the work unhin¬ 
dered. Thus with the view of being helpful to the readers we present the detailed calculation. 


A.l Expressions Related to Spin Oscillation Frequency 

We will start this section with the introduction of the line element, 

ds^ = —f{r)dt^ + -I- r'^dVL^ 

9(r) 


(57) 


Then the Christoffel symbols corresponding to the above metric ansatz are given by: 
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(58) 
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Then we have the following expression for various covariant derivative components of the vierbein 
vectors as 


V.e(o)^ = 

Vte(i)t =Vee(i)e=rv^; = r^/g sin^ 6 

Vee( 2 )r =- 1 ; V 0 e( 2)0 = rsindcosd 

V,^e( 3 )r =-sind; V 0 e( 3 )e =-rcosd. (61) 
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Components of are given by: 


i(0) = v7c^°; = rC/^ = rsin6»C/‘^ 

Va 


The components of Q{a){b) are as follows: 
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Hence components of electric and magnetic fields are: 
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( 68 ) 


If we specify to the 0 = 7r/2 plane, then we readily arrive at the following expressions for the 
components of G(q) as, dO/dr = 0, and consequently 5^ = 0 as. 


G(i) — 0; 


G(3) — 0; 


G(2) = - 


u* 


V~a- 


fVia/f) 
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(69) 


These are the expressions used in the main text while calculating the spin oscillation frequency. 
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A.2 Expressions Related to Flavour Oscillation Frequency 


We will study the flavour oscillation frequency for the general metric ansatz as presented in 
Eq. (57). The conserved energy and angular momentum along with radial momentum has the 
expression: 


= -rrikEk; = rrik^^-, 

9{r) 



TTlkLk 


(70) 


where, E and L, were defined through the relations: 
following function: 


V = 



we arrive at the required expression for f given by: 
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Then the derivatives dt/dr and d(f)/dr can be obtained as, 


(72) 


dt i Ek 

VTgVElv-f 

dip _ Ip _ Lk jj 

dr~ r~ 9 


(73) 

(74) 


Finally the upper component of has the following expression: 


pik)r _ 
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Hence the phase along the geodesic of neutrino, with the assumption that it is massive leads to 
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(76) 


However in the literature sometimes the neutrino though have a mass is taken to travel along 
the null geodesics, as it is extremely relativistic. For that purpose the expressions for f, dt/dr 
and dcp/dr leads to. 



M _ 1 dp _ Lk ly 
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(77) 
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Then the phase along null trajectory can be presented as: 
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B Circular Orbit in Spherically Symmetric Spacetime 

We start with the Lagrangian on the equatorial i.e. 0 = 7r/2 plane, by exploiting the spherical 
symmetry of the problem, such that 


^ 1 .fdty 1 /dry 1 2 fdcj) 


(79) 


Then the energy and the angular momentum are conserved, since the Lagrangian does not involve 
any functions of time or of azimuthal angle (j) [73] . Then we get the equation of motion for t and 
(5 as: 


. E 


since^ (fi) = 0^ 


since— ) = 0 


(80) 

(81) 


From the above expressions it is evident that the four-velocity components are being given by: 

(82) 


’r2 

Thus using the on-mass shell condition: = — 1, we arrive at the equation for r leading to 

the following expression: 


/ 




(83) 


Then along with this relation we can use the expression for (j) leading to the equation for the 
orbit of a massive particle as: 


dr 
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(84) 


Then introducing a new variable r = (1/u), we get, {dr/d(j)) = —{l/v?){du/d4>), which modihes 
the above equation to the form: 
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Differentiating this expression again we obtain the following second order differential equation 
satisfied by the variable u as: 


cPu r g' fl + L^u^\ (g' /'\1 

dcj)^ "" 2g\ ) 2L^u‘^\gf p) 


( 86 ) 


where “prime” denotes derivative with respect to the radial coordinate r. For circular orbit, u 
is fixed, say at u = Uq- Then {cPu/dp) should vanish along with {dr/d(j)). This leads to the 
following equations: 



Then the energy and angular momentum for the circular orbit are obtained as: 


(87) 
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Thus the four-velocity components are being given by: 
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Thus we get the following angular velocity having the expression: 

dt [/* V 2ro 

However for null trajectory i.e. for photons we get equation for orbit as: 

/ dr r'^g f E^ \ 
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the above equation for the orbit can be written introducing the new variable u 
the following form: 
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Then from the condition of circular orbit at r = tq, or equivalently u = Ug, we get the following 
relations between energy and angular momentum such that: 
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For Schwarzschild spacetime, we get the following equation: 2Muq = 2uo(l — 2Muo), leading to 
ro = 3M as the photon circular orbit radius. Then four velocity has the expression: 


[/^ = 




Thus the corresponding angular velocity turns out to be: 


dt ~ m ~ 


(97) 


(98) 


The above expressions for energy and angular momentum, along with the angular velocity have 
been used extensively in the main text. 
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